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Recent analyses that include cosmic microwave background (CMB) anisotropy measurements from 
the Atacama Cosmology Telescope and the South Pole Telescope have hinted at the presence of a 
dark radiation component at more than two standard deviations. However, this result depends sensi- 
tively on the assumption of an HST prior on the Hubble constant, where Hq = 73.8±2.4 km/s/Mpc 
at 68% c.L. Here we repeat this kind of analysis assuming a prior of Ho — 68±2.8 km/s/Mpc at 68% 
C.I., derived from a median statistics (MS) analysis of 537 non-CMB Ho measurements from Huchra's 
compilation. This prior is fully consistent with the value of Ho = 69.7 ± 2.5 km/s/Mpc at 68% c.l. 
obtained from CMB measurements under assumption of the standard ACDM model. We show that 
with the MS Ho prior the evidence for dark radiation is weakened to ~ 1.2 standard deviations. 
Parametrizing the dark radiation component through the effective number of relativistic degrees of 
freedom N^s, we find N^s = 3.98 ± 0.37 at 68% c.l. with the HST prior and N^s = 3.52 ± 0.39 at 
68% c.l. with the MS prior. We also discuss the implications for current limits on neutrino masses 
and on primordial Helium abundances. 
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I. INTRODUCTION 

Recent measurements of the cosmic microwave back- 
ground (CMB) radiation anisotropy made by the Ata- 
cama Cosmology Telescope (ACT) [l| and by the South 
Pole Telescope (SPT) [2[ have provided valuable general 
confirmation of the theoretical predictions of the shape 
of the CMB anisotropy at arcminute angular scales, in 
the diffusion damping regime. 

While the inclusion of these new small-scale data do 
not significantly alter the constraints on parameters of 
the "standard" ACDM cosmological model 0, 0] , com- 
pared to those obtained by using the WMAP satellite 
CMB anisotropy data in conjunction with other cosmo- 
logical measurements Q , they can be used to significantly 
improve the constraints on those new, "beyond-standard- 
modcl" , parameters that mostly affect the physics of the 
CMB anisotropy diffusion damping tail. 

In particular, the recent ACT and SPT data have 
placed new constraints on the number of relativistic de- 
grees of freedom, N^s, that defines the physical energy 
density in relativistic particles today, Prad, given by 



Prad 
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where pj is the energy density of the CMB photons at 
present temperature Tj = 2.726 K (see, e.g., Ref. @]). 
In the standard scenario, assuming three active massless 
neutrino species with standard electroweak interactions, 
the expected value is A'off = 3.046, slightly larger than 
3 because of non-instantaneous neutrino decoupling (see, 
e.g., Ref. 0). 

The new data from ACT and SPT, jointly analysed 
with earlier, large-scale WMAP (and other) data, rule 
out the case of A'off = at high statistical signifi- 
cance. That is, for the first time, CMB anisotropy and 



large-scale structure observations confirm the existence 
of neutrinosQ However, these data also seem to prefer 
a value of Ncs 4, hinting at the presence of an addi- 
tional relativistic component (see Refs. [H, 0): over and 
above the three neutrino species in the standard model 
of particle physics. In particular, some of us, Q, found 
Nas = 4.08lQ'gg at 95% confidence level from such an 
analysis and similar results are presented in Refs. f9l.[l0j. 

These results are significant, since they rather strongly 
suggest that CMB anisotropy data (alone or in conjunc- 
tion with other large-scale cosmological data) indicate 
the presence of some kind of "dark radiation" that is not 
seen in any other cosmological data. They have prompted 
the development of many theoretical models in which A'cff 
is larger than 3, [Tl| . 

Several non standard models related to axions or de- 
caying particles, gravity waves, extra dimensions and 
dark energy [l^ can infact predict a larger value for N^e . 

It is therefore crucial to carefully investigate this re- 
sult, to see if it can be strengthened or weakened by, for 
example, considering a slightly different choice of data. It 
is well known in the literature (see, e.g., Refs. [1,Q) that 
Nes is degenerate with the value of the Hubble constant 
Hq. Assuming a prior on the value of the Hubble constant 
is therefore a key step in the determination of A'cff from 
the data. The prior on the Hubble constant used in most 
recent analyses, labeled HST, is a Gaussian one based on 
the results of Ref. ^ with Hq = 73.8 ± 2.4 km/s/Mpc, 
including systematics. 

While this 3% determination of Hq is certainly im- 
pressive, one might wonder if a slightly different Hubble 



Of course, cosmological big bang nucleosynthesis theory in com- 
bination with the observed light nuclei abundances had pointed 
to this earlier. 



2 



constant prior could change the preference for Ncs > 3. 
There arc several indications that a different Hubble con- 
stant prior could be more appropriate. For instance, a 
number of measurements result in a significantly lower 
value of Hq; e.g., the Ref. [l^ summary value is Hq = 
62.3±4 km/s/Mpc. In addition, a standard analysis, un- 
der the assumption of N^s = 3.046, of CMB data alone is 
able (in a flat universe) to constrain the Hubble constant. 
Recent such analyses yield Hq ^ 70 km/s/Mpc, more 
than one standard deviation away from the HST value. 
For example, the analysis of ACT and WMAP7 data in 
Ref. [3 gives Hq = 69.7 ± 2.5 km/s/Mpc. Clearly, there 
is also observational evidence for a significantly smaller 
value of Hq than the HST estimate. Furthermore, it is 
possible that using a prior with a lower value of Hq could 
result in a iVoff- determined from CMB anisotropy and 
other large-scale data that is consistent with the other 
cosmological N^ff determinations. 

There are very many measurements of Hq, over 
5500 Most recent estimates lie in the interval 60- 
75 km/s/Mpc, with error bars on some individual es- 
timates probably being too small, since these measure- 
ments are mutually inconsistent (this is likely a con- 
sequence of underestimated systematic errors in some 
cases). Clearly, what is needed is a convincing summary 
observational estimate of -ffoEl To date, the best tech- 
nique for deriving such a summary estimate — that does 
not make use of the error bars of the individual mea- 
surements — is the median statistics technique; Ref. [l7l | 
includes a detailed description of this technique. 

The median statistics technique has been used to anal- 
yse a number of cosmological data sets. These include 
Type la supernova apparent magnitude data, to show 
that the current cosmological expansion is accelerating, 
[itI [Tsj : CMB temperature anisotropy data, in one of 
the first analyses to show that these data were consistent 
with flat spatial hypersurfaces, [l^; and, collections of 
measurements of the cosmological clustered mass density, 
in one of the earliest analyses to show that this makes up 
around 25-30% of the current epoch cosmological energy 
budget, [l^l. These successes support the idea that a me- 
dian statistics estimate of the Hubble constant provides 
an accurate summary estimate. 

The median statistics technique has been used thrice 
to analyse Huchra's list (at three different epochs). From 
an analysis of 331 measurements (up to the middle of 
1999), Ref. [l3| found an median statistics summary 
Hq ~ 67 km/s/Mpc; from 461 measurements (up to the 
middle of 2003), and from 553 measurements (up to early 
2011), Refs. [Ill, [22| both found a median statistics sum- 
mary Hq = 68 km/s/Mpc. While the estimated statisti- 
cal error bar (given by the scatter in the central Hq val- 
ues) has decreased as the sample size has increased, the 



^ See cfa-www. harvard. cdu/r^huchra/. 

^ And not just for the case at hand, but for many different cos- 
mological parameter analyses, see, e.g., Refs. [l^ . 



larger (and dominant) systematic error bar (estimated 
from the scatter in the summary values of Hq determined 
by different techniques) has changed much less. 

For our analyses here we estimate Hq using the method 
of Ref. [l^l but now excluding from the Huchra list 
of 553 measurements the 16 Hq measurements derived 
from CMB data assuming A^cft = 3.046. We exclude 
these 16 CMB measurements as we want an external 
and independent prior on Hq to use in our analysis 
of the latest CMB datasets. From a median statis- 
tics analysis of the 537 non-CMB measurements we find 
Hq ~ 68 ±2.8 km/s/Mpc (one standard deviation error), 
identical to that found in Ref. [l^l from an analysis of the 
553 measurements. In what follows we refer to the Gaus- 
sian prior based on this value as the median statistics 
(MS) Hq prior. 

Our goal here is to discuss the implications of assum- 
ing the MS prior for Hq, instead of the usual HST prior, 
for current CMB and large-scale structure parameter in- 
ference. We focus much of our attention on the value 
of TVoff and the evidence for dark radiation, but we also 
consider how the MS prior changes the estimated value 
of other parameters, including the dark energy equation 
of state parameter w and the spectral index of primordial 
fluctuations n^. 

Our paper is organized as follows. In the next Section 
we briefly summarize the data analysis method we use. 
In Sec. HI we present our results. We conclude in Sec. 
IV. 



II. ANALYSIS METHOD 

Our analysis is based on a modified version of the 
public COSMOMC [13 Monte Carlo Markov Chain 
code. We include the following CMB data: WMAP7 [1, 
ACBAR [13, ACT [H, and SPT [|, including measure- 
ments up to maximum multipole number lma.K = 3000. 
As in Ref. Q we include galaxy clustering data from the 
SDSS-DR7 luminous red galaxy sample [l^]. Also, as dis- 
cussed in the Introduction, we choose two different priors 
on the Hubble constant: the median statistics (MS) prior 
of Hq = 68 ± 2.8 km/s/Mpc as well as, for comparison, 
the HST prior [TJ] used in previous analyses. 

In the basic analysis we sample the usual seven- 
dimensional set of cosmological parameters, adopting flat 
priors on them: the baryon and cold dark matter den- 
sities U\^h^ and Vt^h^ , the ratio of the sound horizon 
to the angular diameter distance at decoupling 6, the 
optical depth to reionization t, the scalar spectral in- 
dex ng , the overall normalization of the spectrum Ag at 
k = 0.002 Mpc^^, the effective number of relativistic de- 
grees of freedom iVoff. 

Our analysis is very similar to the one presented in 
Ref. [1], with three changes: (i) we consider two differ- 
ent Hq priors; (ii) we consider an extended case where 
we assume massive neutrinos, we enlarge our parame- 
ter space varying the total mass of neutrinos ^ rrii, ; (iii) 
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we allow the Helium abundance Yp to vary consistently 
with standard BBN following Ref. 01- This means that 
each theoretical CMB angular spectrum is computed as- 
suming a value for Yp derived by BBN nucleosynthesis 
from the input values of ilt,h^ and iVeff of the theoretical 
model considered. The small uncertainty on Yp derived 
from the experimental errors on the neutron half-life pro- 
duces negligible changes in the CMB angular spectra so 
we ignore it. In a latter case we also vary Yp as a free 
parameter. 

In addition, where indicated, we also present con- 
straints on the dark energy equation of state parame- 
ter w (the ratio of the pressure to energy density of the 
dark energy fluid), assumed to be redshift independent, 
although the corresponding dark energy density is time 
dependent We consider massless neutrinos, adiabatic 
initial conditions, and a spatially-flat universe. 

Following Ref. Q we account for foreground contribu- 
tions by marginalizing over three additional amplitudes: 
the Sunyaev-Zeldovich effect amplitude Asz, the ampli- 
tude of clustered point sources Ac, and the amplitude of 
Poisson-distributed point sources Ap. 



III. RESULTS 

A. Neutrinos 

As in Ref. Q, we compute the likelihood function in 
the seven-dimensional (or eight when massive neutrinos 
are considered) cosmological parameter space described 
above, and multiply it by the prior probability distri- 
bution functions to derive the seven-dimensional poste- 
rior probability density distribution function. Marginal- 
izing this over all but one of the cosmological parameters 
gives the one-dimensional posterior probability distribu- 
tion function for the parameter of interest. This one- 
dimensional distribution function is used to determine 
the most likely value of the parameter, as well as limits 
on it. These are listed in Table I, for three different Hub- 
ble consta nt p riors: a flat one (no prior); the Gaussian 
HST one, and the Gaussian MS one. Marginalizing 
over only five of the cosmological parameters, we derive 
the two-dimensional posterior probability density distri- 
bution function P{Ho, Ncs). This is used to derive the 



constraint contours in the two-dimensional NoS-Hq pa- 
rameter space shown in Fig. 1, for the two Gaussian Hq 
priors. 



* This is the widely-used XCDM parametrization of dark energy. 
It is not a complete parametrization, as it cannot describe the 
evolution of spatial inhomogeneities, nor is it an accurate approx- 
imation of more physically motivated time-varying dark energy 
models, |26^ . It is preferable to use a consistent and physi cally 
motivated dark energy model, e.g., that proposed in Refs. l27ll . 
for such an analysis, but this is a much more involved under- 
taking, so instead we patch up the XCDM parametrization by 
assuming that the acoustic spatial inhomogeneities travel at the 
speed of light. This extended XCDM parametrization should 
provide reasonable (qualitative) indications of what might be 
expected in a consistent, physically- motivated model of time- 
varying dark energy. 




FIG. 1. Constraints in the NcS-Hq plane. Elliptical two- 
dimensional posterior probability distribution function con- 
tours show the 68% and 95% c.l. limits. Red contours and 
regions (closer to the upper right corner) assume the HST 
prior with Ho = 73.4 ± 2.4 km/s/Mpc, while blue contours 
and regions (closer to the lower left corner) are obtained us- 
ing the median statistics prior with Hq = 68±2.8 km/s/Mpc. 
The dotted black vertical line corresponds to A'^eff = 3.046. 

Table I and Fig. [T] show that the Hq prior plays a cru- 
cial role in determining constraints on Ngg from the data. 
With the HST Hq prior we find a central Nes value that 
is 2.5(7 larger than 3.046, while the median statistics prior 
results in an A'cff that is consistent with 3.046 (being only 
1.2(7 larger). 

The HST prior is therefore at least partially responsi- 
ble for the current indication for dark radiation. How- 
ever, as we can see from the central values of Hq and 
A'cff obtained when a flat prior on Hq is assumed, the 
CMB anisotropy and large-scale structure data consid- 
ered here prefers a larger value of A'eff (being 1.9a larger 
than 3.046) and a somewhat larger value of Hq. This is 
clear also from the Xmin values of the best fit that are 
higher when the median statistics Hq prior is assumed, 
compared to the case of the HST prior (see the last line 
of Table I). 

The Hq prior is crucial also in the determination of the 
^ rrLiy limits if we instead limit ourselves to the case of 3, 
standard, massive neutrinos. In Table HI columns 3 and 
5, we quote the cosmological parameters and the upper 
limits on ^ rrii, in case of the HST and of the MS prior. 
As we can see, the upper limit on ^ rrii, is considerably 
weaker when the MS prior is considered, with the 95% 
c.l. upper limit moving from ^TOj^ < 0.36 eV in the 
case of the HST prior to < 0.60 eV in the case 

of the MS prior. This can be clearly explained by the 
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Parameters 


No Prior 


HST Prior 


MS Prior 






73.8 ± 2.4 km/s/Mpc 


68 ±2.8 krn/s/Mpc 




0.02258 ± 0.00050 


0.02248 ± 0.00039 


0.02210 ± 0.00037 


0.02211 ± 0.00040 


0.02188 ± 0.00036 




0.134 ±0.010 


0.1317 ± 0.0080 


0.1142 ± 0.0029 


0.1256 ± 0.0080 


0.1181 ± 0.0032 


9 


1.0395 ±0.0016 


1 0397 ± 0016 


1.0415 ± 0.0014 


1 0400 ± 0017 


1 0409 ± 0014 


T 


0.085 ±0.014 


084 ± 013 


083 ± 013 


080 ± 013 


081 ± 014 


Us 


0.984 ± 0.017 


0.979 ± 0.012 


0.9659 ± 0.0091 


0.964 ± 0.012 


0.9536 ± 0.0090 


iVcff 


4.14 ±0.57 


3 98 ± 37 


3.046 


3 52 ± 39 


3.046 


Em.[eV] 


0.0 


0.0 


< 0.36 


0.0 


< 0.60 


Ho [km/s/Mpc] 


75.2 ±3.6 


74.2 ± 2.0 


69.3 ± 1.4 


70.9 ± 2.1 


66 8 ± 1 8 

\J\J tKj _]_ _L ■ 


log{10^°As) 


3.183 ±0.043 


3.191 ±0.035 


3.205 ± 0.034 


3.219 ±0.036 


3.226 ±0.034 




0.277 ±0.019 


0.280 ±0.016 


0.284 ±0.017 


0.294 ±0.017 


0.315 ±0.024 




0.882 ± 0.033 


0.876 ± 0.028 


0.782 ± 0.032 


0.857 ± 0.028 


0.757 ± 0.043 


^SZ 


< 1.4 


< 1.3 


< 0.97 


< 1.1 


< 0.96 




< 14.5 


< 14.7 


< 12.8 


< 14.1 


< 13.1 




< 24.9 


< 25.5 


< 26.6 


< 26.1 


< 26.6 


Xmin 


7593.4 


7593.2 


7592.0 


7594.8 


7595.1 



TABLE I. Cosmological parameter values and 68% confidence level errors assuming A'cff relativistic neutrinos or A'^eff = 3.046 
massive neutrinos. 95% c.l. upper bounds are listed for the sum of neutrino masses and foregrounds parameters. We also list 
the derived Hubble constant, the non-relativistic matter density parameter fim = ilc + fit, and erg, the amplitude of density 
inhomogeneities averaged over spheres of radius 8/i~^ Mpc, where h is the Hubble constant in units of 100 km/s/Mpc. 



CMB degeneracy between and ^ nii, as illustrate in 
Fig. [21 Namely, lower values of the Hubble parameter 
are in better agreement with current CMB data when 
^ rrij/ is increased. Dataset preferring higher values for 
i?o will therefore provide stronger constraints on ^ rriv 
when combined with the CMB data. 

Beside the NcS-Hq degeneracy, it is interesting to note 
that there also is a degeneracy between Ngg and Ug. 
When the HST prior is assumed, Ug is l.Sa below 1, 
while for the median statistics case it is 3(t below unity. 

In Fig. Owe show the contours in the two-dimensional 
^m~crs parameter space, for the two Gaussian Hq priors. 
Here erg is the amplitude of density inhomogeneities av- 
eraged over spheres of radius 8h~^ Mpc. In this figure we 
also show the fit to the central value and the two standard 
deviation limits of the constraint from the normalization 
of the galaxy cluster mass function from Ref. [1^, i.e., 
as = (0.25/17m)° ''^[0.813± 0.013 ±0.024]. Here the first 
error bar represents the statistical, and the second the 
systematic, error (see their Sec. 10). We derive the 2a 
cluster constraints shown in Fig. 2 by adding these errors 
in quadrature and then doubling. 

From Fig. [3] we see that both Hq priors give results that 
are not far off from what the measured normalization of 
the cluster mass function demands. Qualitatively, the 
HST Hq prior is more consistent with the cluster data if 
flm ~ 0.25, near the low end of current indications, see, 
e.g., Ref. [201, while the median statistics case prefers a 
larger J^ni ^ 0.27, more consistent with current measure- 
ments, see, e.g., Ref. po| . 




'0 0.2 0.4 0.6 0.8 

Zm [eV] 



FIG. 2. Constraints in the m^-Ho plane. Elliptical two- 
dimensional posterior probability distribution function con- 
tours show the 68% and 95% c.l. limits. Red contours and 
regions (closer to the upper left corner) assume the HST prior 
with Ho = 73.4 ± 2.4 km/s/Mpc, while blue contours and re- 
gions (closer to the lower right corner) are obtained using the 
median statistics prior with Ho = 68 ± 2.8 km/s/Mpc. 
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FIG. 3. Constraints in the Q„^-as plane. Elliptical two- 
dimensional posterior probability density function contours 
show the 68% and 95% confidence level limits. Red contours 
(closer to the upper left corner) assume the HST prior with 
Ho ~ 73.8±2.4 km/s/Mpc; blue contours (closer to the lower 
right corner) are obtained with the median statistics prior 
where Ho = 68 ± 2.8 km/s/Mpc. The green region (in the 
lower left corner) demarcates the central value and 2a limits 
from the cluster mass function normalization data, [2^ . 



B. Helium mass abundance 

One assumption made in the previous paragraph is 
that the Helium abundance is varied consistently with 
BBN. Current CMB data produce only weak constraints 
on this quantity and allowing Yp to vary freely would 
make the standard case of iVoff — 3.046 in better agree- 
ment with data due to an anti-correlation between iVoff 
and Yp in CMB data (see, for example, the discussion 
in [HI). In order to check the impact of the Hq priors 
in this case, we have performed two analysis varying the 
Helium abundance Yp and N^ff . The results are reported 
in Table EL 

As we can see, when Yp is allowed to vary, the stan- 
dard case of Naff is more consistent with current data 
in both cases. In the case of the MS prior we have 
Ncff = 2.75 ± 0.46 that is perfectly consistent with the 
expectations of the standard scenario. However the value 
obtained for the Helium abundance is probably too high 
in the case of the MS prior: Yp = 0.334 ± 0.033 that is 
about two standard deviations away from the conserva- 
tive experimental bound of 1^ < 0.2631 obtained from 
an analysis of direct measurements in [30l |. 

The larger helium abundance obtained in the case of 
the MS prior respect to the HST prior can be clearly seen 
from the direction of the degeneracies in the 2D contours 
plots in Figure H] Namely, a lower iVoff prefers an higher 
Yp and a lower prior for Hq shifts the constraints towards 
lower N^ff and higher values for Yp. 



Parameters 


HST Prior 


MS Prior 




0.02274 ± 0.00042 


0.02246 ± 0.00043 




0.1246 ±0.0091 


0.1138 ± 0.0085 


e 


1.0429 ± 0.0027 


1.0454 ± 0.0029 


T 


0.087 ±0.014 


0.085 ±0.014 


ris 


0.986 ±0.013 


0.972 ±0.013 




3.52 ±0.48 


2.75 ± 0.46 


i/o [km/s/Mpc] 


72.7 ±2.2 


68.2 ±2.3 


Yp 


0.310 ±0.034 


0.334 ± 0.033 




3.175±0.037 


3.197 ±0.036 




0.279 ±0.015 


0.293 ±0.016 




0.872 ± 0.029 


0.847 ± 0.029 


Asz 


< 1.7 


< 1.6 




< 15.4 


< 15.3 




< 23.1 


< 23.4 


Xmin 


7592.0 


7590.4 



TABLE II. Cosmological parameter values derived assuming 
a varying Yp. Errors are at 68% c.l. while upper bounds at 
95% c.l. are reported for foregrounds parameters. 



C. XLCDM 

The standard ACDM cosmological model has some 
conceptual problems that are partially alleviated in some 
models in which the dark energy density varies slowly in 
time (and so weakly in space), [l^l- Furthermore, ob- 
servational constraints on cosmological parameters are 
model dependent, i.e., the observational estimate of a 
cosmological parameter, e.g., N^ff, depends on the cos- 
mological model used to analyse the data. It is therefore 
of interest to examine the observational cosmological con- 
straints on A^off in a cosmological model in which the dark 
energy density varies in time, such as that of Ref. [l^] . 
This is a somewhat challenging task that we will leave 
for future work. However, to get an indication of what 
could be expected from such an analysis, we determine 
the observational constraints on A"off in a cosmological 
model in which the time-evolving dark energy density 
is parametrized by the XCDM parametrization (made 
complete by assuming that the acoustic spatial inhomo- 
geneities propagate at the speed of light) described above. 
Table Hill shows the observational constraints derived un- 
der these assumptions. 

From Table Hill we see that the MS prior changes the 
best fit w in the standard case with A"off = 3.046 to 
w ^ —0.9, with w = —1 off by one standard deviation. 
When both if and A'ofr arc allowed to vary freely the geo- 
metrical degeneracy with Hq makes the HST and MS Hq 
priors much less effective. In this case the evidence for 
dark radiation is again significant: for the MS Hq prior 
case we find A^eff = 4.16 ± 0.53, and a dark energy equa- 
tion of state parameter w = —0.76 ± 0.10, i.e., excluding 
a cosmological constant at more than two standard de- 
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Parameters 


HST Prior 


MS Prior 




0.02200 ± 0.00040 


0.02290 ± 0.0054 


0.02206 ± 0.00040 


0.02279 ± 0.00053 




0.1162 ± 0.0039 


0.1347 ±0.0085 


0.1141 ±0.0040 


0.1291 ± 0.0084 


e 


1.0414 ±0.0015 


1.0396 ±0.0016 


1.0414 ±0.0015 


1.0400 ± 0.0016 


r 


0.080 ±0.013 


0.089 ±0.015 


0.081 ±0.013 


0.089 ±0.015 


Us 


0.956 ±0.010 


0.997 ±0.019 


0.959 ±0.011 


0.993 ±0.019 


iVcff 


3.046 


4.42 ± 0.54 


3.046 


4.16 ±0.53 


Ho[km/s/Mpc] 


72.1 ± 2.4 


72.8 ± 2.3 


66.7 ±2.6 


68.0 ± 2.4 


w 


-1.09 ±0.10 


-0.86 ±0.11 


-0.90 ±0.10 


-0.76 ±0.10 


logiW^'^As) 


3.223 ± 0.039 


3.150 ±0.050 


3.210 ±0.041 


3.149 ±0.051 




0.267 ±0.018 


0.298 ± 0.022 


0.307 ± 0.226 


0.329 ± 0.025 


0"8 


0.856 ± 0.044 


0.831 ±0.047 


0.790 ± 0.046 


0.775 ± 0.047 


Asz 


< 0.94 


< 1.5 


< 0.95 


< 1.4 




< 13.0 


< 15.0 


< 13.0 


< 14.9 


Ap[tiK^] 


< 27.0 


< 23.9 


< 26.7 


< 24.7 


Xmin 


7598.1 


7592.7 


7595.1 


7592.1 



TABLE III. Cosmological parameter values derived assuming the XCDM parametrization of time-evolving dark energy. Errors 
are at 68% c.l. while upper bounds at 95% c.l. are reported for foregrounds parameters. 



Parameters 


HST Prior ± SNela 


MS Prior + SNela 




0.02203 ± 0.00038 


0.02260 ± 0.00046 


0.02190 ± 0.00038 


0.02230 ± 0.00046 




0.1156 ±0.0037 


0.1317 ±0.0079 


0.1157 ±0.0038 


0.1249 ±0.0077 


e 


1.0414 ±0.0015 


1.0400 ± 0.0016 


1.0411 ±0.0015 


1.0401 ±0.0016 


T 


0.081 ±0.013 


0.086 ±0.014 


0.080 ±0.013 


0.083 ±0.014 


ris 


0.957 ±0.010 


0.985 ± 0.015 


0.956 ±0.010 


0.972 ±0.016 


N,« 


3.046 


4.08 ± 0.43 


3.046 


3.63 ±0.42 


ffo[km/s/Mpc] 


71.0 ± 1.6 


74.0 ± 2.0 


68.8 ± 1.6 


70.6 ± 2.1 


w 


-1.050 ±0.069 


-0.967 ± 0.075 


-0.989 ±0.070 


-0.946 ± 0.076 


/o3(10"A,) 


3.222 ±0.038 


3.178 ±0.043 


3.221 ±0.038 


3.198 ±0.044 




0.273 ±0.014 


0.282 ±0.015 


0.291 ±0.015 


0.295 ±0.016 




0.843 ± 0.035 


0.863 ± 0.038 


0.823 ± 0.036 


0.836 ± 0.038 


Asz 


< 0.94 


< 1.3 


< 0.94 


< 1.2 


Ac[^K^] 


13.0 


14.8 


< 13.1 


< 14.0 


Ap\^iK^] 


27.0 


24.8 


< 27.0 


< 26.0 


Xmin 


8128.4 


8124.0 


8126.2 


8125.6 



TABLE IV. Similar constraints as in Table Hill but now also including the SNela data in the analysis. 



viations. A scale-invariant HPYZ primordial spectrum 
with = 1 is fully consistent with both priors. While 
some of these values indicate significant tensions with the 
standard ACDM model, it is important to keep in mind 
the strong degeneracies between N^e, Hq and w, as well 
as the fact that the XCDM parametrization used in the 
analysis has been arbitrarily completed to allow for an 
accounting of the evolution of density inhomogeneities. 

In order to try to break these degeneracies, and derive 
more reliable constraints on the parameters, we perform a 
new analysis that also include the SDSS supernova Type 
la (SNela) apparent magnitude data, [l^. From Table 



IIVI we see that the inclusion of the SNela data bring the 
results back to the previous dichotomy: the HST prior 
clearly shows a preference for N^s > 3.046 while the MS 
prior results in a value of A'ofr that is in much better 
agreement with the standard scenario. The constraints 
on the equation of state are w = — 0.967 ± 0.075 for the 
HST prior and w = -0.946 ±0.076 for the MS prior. The 
HPYZ spectrum with = 1 is again in tension with the 
observations for the MS Hq prior at a little less than two 
standard deviations. 
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IV. CONCLUSIONS 

A "standard" cosmological model is starting to fall in 
place. Interestingly, recent data have provided some in- 
dication for an unexpected new "dark radiation" com- 
ponent. In this brief paper we have again emphasized 
the important role played by the HST Hq prior in estab- 
lishing the statistical evidence for the existence of this 
dark radiation. We have also shown that with a new me- 
dian statistics Hq prior derived from 537 non-CMB Hq 
measurements, there is no significant evidence for A^off > 
3.046, consistent with the indications from other cosmo- 
logical data. And it is probably not unreasonable to be- 
lieve that the converse might also be true: with other cos- 
mological data not inconsistent with A^eff = 3.046, con- 
sistency of the smaller-scale CMB anisotropy data with 
the predictions of the ACDM model apparently demands 
Ho - 68 km/s/Mpc. 

We emphasize, however, that when the same data 
are analysed in the context of the somewhat arbitrarily- 
completed XCDM dark energy parametrization, they 
prefer N^s > 3.046. It would be useful to see if this 
remains the case if, instead of XCDM, a complete and 
consistent model, such as (pGDM, is used in the analysis. 

We have shown that the HST Hq prior is, at least par- 
tially, responsible for the evidence supporting the exis- 
tence of a new dark radiation component. However, fu- 
ture CMB anisotropy and galaxy clustering data, as well 
as a definitive determination of Hq, will be needed to 
fully resolve this issue. 

In particular, the future data expected from the Planck 
satellite should be able to constrain independently the 
values of Hq and A^eff, clarifying if the current tension 
between the HST and CMB constraints on Hq is due to 
a dark radiation component or systematics in the data. 



FIG. 4. Constraints in the Yp - N^ff plane (top) and Yp - 
Ho (bottom). Elliptical two-dimensional posterior probabil- 
ity distribution function contours show the 68% and 95% c.l. 
limits. Red contours and regions (closer to the upper left cor- 
ner) assume the HST prior with Hq = 73.4 ± 2.4 km/s/Mpc, 
while blue contours and regions (closer to the lower right 
corner) are obtained using the median statistics prior with 
Ho = 68 ± 2.8 km/s/Mpc. 
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